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Correlation functions of the Lieb-Liniger gas and the LeClair-Mussardo formula 


Garry Goldstein and Natan Andrei 
Department of Physics, Rutgers University and 
Piscataway, New Jersey 08854 

In this letter we derive formulas for multi point correlation functions, in the thermodynamic 
limit, for the Lieb Liniger gas taken with respect to arbitrary eigenstates. These results apply for 
the ground state, thermal states and GGE states. We obtain these correlation functions as a series 
of multiple integrals of progressively higher dimensions. These integrals converge rapidly for short 
distance correlation functions and low densities of particles. The series derived matches exactly the 
LeClair Mussardo formula for correlation functions of relativistic integrable models. 


Introduction. The field of integrable models has 
reached a certain level of maturity. For many integrable 
models in 1-D the exact eigenstates and eigen energies 
are know for arbitrary size systems and their form in the 
thermodynamic limit is well known through the string 
hypothesis The thermodynamics of many inte¬ 

grable models is well understood through the thermo¬ 
dynamic Bethe ansatz which is also based on the string 
hypothesis EHi- One of the key remaining challenges is 
the calculation of correlation functions for local observ¬ 
ables 0. This task is very important as the theoretical 
description of a many body system is usually given in 
terms of correlation functions of local observables. In 
particular the results of any experiments that may be 
carried out on the system can be computed in terms of 
sufficiently complicated multi point correlation functions, 
G{xi ■ ■ ■ Xn) = (Oi(xi) • • • On{xn)), with the expectation 
value is taken with respect to the initial state - typically 
the ground state or a thermal ensemble (which is known 
to be equivalent to an eigenstate calculation i). Many 
correlation functions are directly measurable. For exam¬ 
ple, in the cold atoms context, correlation functions are 
measurable through time of flight interferometry (Ml. 
The calculation of these correlation functions has been a 
major challenge over the wars and many approximation 
methods were developed (2, [3 ■ Even when the system is 
described by an integrable Hamiltonian with all its eigen¬ 
states given by the Bethe Ansatz d or by muivalently 
by the Quantum Inverse Scattering method (2|, the cor¬ 
relation functions are known for only few cases and are 
implicitly expressed in terms of complicated infinite de¬ 
terminants. 

The study of correlation functions for arbitrary states, 
rather than with respect to the ground state, has recently 
been motivated by the study of quench dynamics in inte¬ 
grable systems. In some cases - when no bound states oc¬ 
cur and when the initial state is translationally invariant 
or close to it - the system is described in the long time 
limit by a generalized Gibbs ensemble (GGE) [isl - l^ . 
Eurthermore the equilibrium state which corresponds to 
a GGE ensemble is equivalent to a specially chosen eigen¬ 
state of the integrable model [23 ■ By judicious choice 
of the initial state virtually any eigenstate may be ob¬ 
tained at long times. It is of great theoretical interest 
to compute correlation functions for such eigenstates (or 


equivalently understand quenches at long times). 

The model we will study the Lieb-Liniger Hamil¬ 
tonian describes bosons moving on the continuous 
line and interacting via a short range potential 
[23 ■ Imposing periodic boundary condition with pe¬ 
riodicity L the Hamiltonian is given by, Hrr = 

I-i ^/2 i^) {x) + c (cc) b (a:))^|, with b^x) 

being the creation operator of the bosons at point x. The 
IV-bosons eigenstates of the model, \ki...kN), labeled by 
rapidities {ki.-.k^} = {fc}, are explicitly given by: 


I^Ll2dxN.:J^[%dxiX 

X ni<j (ki -kj)ll n b"^ i^i) |0) 


( 1 ) 


with energy E = ^kf and momentum P = '^ki. The 
scattering factor (k) = incorporates the 

S'-matrix of the Lieb-Liniger model, 

The S'-matrix describes collisions between two particles 
in the model and the integrability of the model hinges 
on the fact that products multi-particle collisions can 
be consistently described in terms of 2-particles col¬ 
lisions [IM3- " The correlation functions in the Lieb- 
Liniger model take the form ({A:}| Oi{xi) ■ ■ ■ On(xn) |{fc}) 
where the states |{fc}) will be specified by the den¬ 
sity Pp{k) of the Bethe-Ansatz momenta k and the 
density of the holes ph{k). Both densities are ob¬ 
tained from solving the Bethe-Ansatz equations that 
follow from imposing periodic boundary conditions, 
exp(iLfcj) Yik^j ^LL (kk — kj) = I. In particular all cor¬ 
relation functions will depend explicitly on the occupa¬ 
tion probability, / (fc) = with ptik) = Pp{k)+ph{k) 

being the total quasiparticle density. Such states |{fc}) 
capture, for appropriately chosen states both thermal av¬ 
erages or their GGE generalizations which are the long 
time limit of the Lieb-Liniger system quenched from some 
given initial state. Thus the ability to compute expecta¬ 
tion values is very useful in understanding the not only 
the thermodynamics of this model (and other integrable 
models) but also the non equilibrium evolution dynamics 
of integrable many body quantum systems. 

In this letter we present a series expansion formula 
for arbitrary multipoint correlation functions for arbi¬ 
trary eigenstates of the Lieb Liniger gas. This expansion 
is given in terms of a series of finite dimensional inte- 
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grals where all the information about the eigenstate be¬ 
ing given in terms of the occupation probability for the 
state. The series converges efficiently for a small density 
of particles (or equivalently a large coupling constant) 
and for short distances between the points of the correla¬ 
tion function. The series presented is identical in form to 
the LeClair Mussardo formula for correlation func¬ 
tions of relativistically invariant integrable models. More 
precisely we will prove that for a generic local operator O 
(which can be a multipoint functions) for the Lieb liniger 
gas we have that: 


1 p ^ jp, 

{{k}\o\{k}) = ■•*") 

n ■ j=l 

( 2 ) 


Here F^„ (fci,..fc„) is appropriately chosen and depends 
only on fci...fc„ (in particular there are no dressing equa¬ 
tions (29l - [3^ ). We will derive explicit formulas for 
F^{ki,..kn) for short distance expansions of the field 
field correlation function, for the density-density corre¬ 
lation function and for the emptiness probability all of 
which are directly measurable in cold atoms experiments. 

Algebraic Bethe ansatz. The correlation functions are 
conveniently expressed in the equivalent language of the 
Algebraic Bethe Ansatz (ABA) as it applies to the Lieb 
Liniger model. The main object used in the ABA is the 
so called transfer matrix [^: 


T{k) 


A{k) B{k) \ 
C {k) D {k) ) 


( 3 ) 


with Aik), B {k), C (k) and D (k) are operators on the 
space of the bosons, e.g. A{k) = A {k, {6 (x)} , {b^ iv)}) 
etc. Also C (k) = B^ (fc). The vacuum eigenvalues of the 
operators A{k), D {k) and C (fc) are given by: C {k) |0) = 
0, A (k) |0) = a (k) |0) , D (k) = d (k) |0), with a (k) = 
exp(— d{k) = exp(i|^fc). It is possible to show 
that the state B (fci) ...B (/cat) |0) is an eigenstate of the 
Lieb-Liniger Hamiltonian if the rapidities {ki} satisfy the 
Bethe Ansatz equations. It is convenient to normalize our 
states and write: 


the velocity distribution and is directly measurable in 
experiment. It is known that the correlation function for 
b^ (x) b (y) with respect to a state |{fc}) with occupation 
density f (k) = is given by [J: 

{{k}\ 6t (x) b jy) life}) ^ (0| det (1 + ^Wt) |0),^o 
{{k} I {k}) det (I - ^Kt) 

( 5 ) 

with the determinant extending over all the rapidities km 
entering the state |{/c}) and where the operator valued 
matrix is given by: 

Wt {k, q) [t {k, q) + t {q, k) exp (z (y - x) {q - k) + 

+^^2 (?) - (?) + {k) - ^A2 {k)) - 

- [t {k, q) exp (z (y -x)[q-k) + (fc) -f 
(g)) + t (q, k) exp ($^, (y) -k (A:))] 
X exp {k) - V'Ai {k) - (fc) + 

+■0212 (g) + 0^3 (g) - 0 Di (g) - ^^2 (g)) + 

-I- acexp (—z (y — x) fc -I- (fc) — $^12 (^) + 

+ 0 A 3 (g) + 0 a 2 (g) - 0 rii (g) - (g))] 

( 6 ) 


Here f {k,q) = ^ , g{k,q) = h{k,q) = 


and t{k,q) = ■ The defining relations for 

these fields are given by: = Pa^ + Qa^ , = 

Pok + Qd^ and V'A, = PAk + dAu , = PDk + qo^ ■ 

The various expectation values are fyj |0) = Pa |0) = 
(0| qa = (6| Qa = 0, (0 I 0 ) = 1. Here |0) is not re¬ 
lated to any states of the boson system but is a ficti¬ 
tious auxiliary ground state. The only nonzero commuta¬ 
tion relations are given by: = Sik In {h (y, fc)), 

[PDi,QDk] = Sik\n{h{k,q)), [pAi,qAk] = 5ikln{h {q,k)) 
and [pDi,qDk] = In (h (A:, y)). Furthermore the de¬ 
terminant Wt {k, q) in Eq. ([6]) is well defined since 
(A;), $6 (y)] = [0a (A:), 0h (y)] = [$a (A:), 0h (g)] = 0. 
Using the Taylor series formula for Fredholm determi¬ 
nants [U it is possible to obtain a series solution for the 
determinant: 


\kT..kN) = n |0) 

j<k k 

where with f{k,q) = . With this normalization 

it is possible to show that the wavefunction of the 
state \ki,...kj^) is given by Eq. ([T]) above. Eurthermore 
the states are normalized such that: 

{ki...kN I ki....kN) = det {Mjk) (4) 

with Mjfc = Sjk (l + Ei=i “ c^+W-kkf^ 

see (21 • 

Field field correlation functions. As an example of the 
general formalism we proceed now to calculate the field 
field correlation functions. Its Eourier transform yields 


^ (0| det (1 -f ^Wt {k, y)) |0) = 

= 0 ^ / fr n / ikj) det (0| (Wt {k„ kk)) |0) 

( 7 ) 

Here ( Wt (kj, kk) ] is the nxn matrix whose j, A;’th 

V /nxn 

entry is Wt {kj,kk), where WT(A:,y) = Simi¬ 

larly we have that 

(1 - ^Kt) = (8) 

= Er=o / $ n / (tj) (il-rlfe.) 

Hence we obtain a LeClair Mussardo like formula for 
the correlation function by Taylor expanding the ratio 
of these determinants, 
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(6|det {1 + ^Wt)\0) 
det (1 - ^Kt) 




n=0 


pb^{x)b{y) 


{ki, ..kn) 


where the {ki, ..kn) are given by 

pbH.Wv) ^ ^det (0| (Wr ikj,k,))^^^ |0) 

pbHcoMy) ^ + ( 6 | (WTikj,kj))^^^\Q)-(^kTikk,kk))^^^ 


etc and where (Kt {kj, kk)\ is the n x n matrix whose j, fc’th entry is Kt (kj, kk) = calculate 

the first few terms in the expansion for the field field correlation function, it is given by: 

(fet (x)5(y))2 = j dk^^ exp [-i{y-x) k]dk + J dk^^ J dq^^‘^exp[-i {y - x) k] + 

j ~ 

+ * J j exp [-i {y - a;) q] (^-2i {x - y) f~^ {k, q) + (g, k) + (fc, + .... 

( 9 ) 


From this expression we see that the LeClair Mussardo 
formula is a short distance expansion (in particular the 
term ^ x — y diverges at large distances). This is a 
general feature of the expansion. This expression is valid 
for X > y, while for x < y one needs to replace the 
expression by its complex conjugate. In the supplement 
we give some examples of how to use these formulas for 
various initial states. 

Density density correlation functions. We now pro¬ 
ceed to consider density density correlation functions, 
this would be another completely explicit example. Ac¬ 
tually it is easier to consider the generating function 
of density density correlation exp (aQxy), where Qxy = 
f^bHz)b{z)dz. It is known that the correlation func¬ 
tion for exp(Q:Qa;y) with respect to a state |{fc}) with 
occupation density / (fc) = is given by 0 : 

({fell exp jaQxy) |{fc}) _ ( 6 | det (l -f ^Vr) |0) 

{{k}\{k}) det{l-^Kr) ^ ^ 

Here the we have the ratio of two Fredholm determinants. 

(k, q) = (^k-qf^+c^ f ( 9 ): and 

Vt (fc, q) =IS^ [t (fc, q) + t (g, k) exp {-i {y - x) {k - q)) x 
c 

X exp {ipi iq) - ipi (/c)} -f e“ exp {(^3 (fc) -f 

+<P 4 (q)} X {t [q, k) + t {k, q) X 

exp {-i {y -x){k- q)) exp {vJ 2 {k) - (g)})] 

( 11 ) 


The various fields are defined as ipi (fc) = qi (k) -I- 
Pi {k), i = 1, 2,3,4. The commutation relations for these 
fields are given by [pa (k) ,pb (q)] = [qa (k), qb (q)] = 0 
and [pa {k), qb (q)] = Ha,b{k,q). With Ha,b {k,q) = 
(In [h [k, q))) A + (In (h [q, k))) A^, 


/ -1 0 0 -1 \ 
0-11 0 I 

1 0-11 
Vo -1 1 - 1 / 


The vacuum states |0) and (0| are defined as pa {k) |0) = 
(0| qa (k) = 0, (0 I 0 ) = 1. Furthermore this determinant 
is well defined since [ipi{k),(pj (q)] = 0, i,j = 1,2, 3,4. 
Proceeding with a Taylor expansion much like in the field 
field case we obtain that the density density correlation 
function is given by: 
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(pW?(!/)) =(exp(afe,)) = 

= i (6| {2e“2t („, k) t (k. «) exp (-j („ - x) (k - ,)) X 

X exp ((/ 5 i {k) - ipi (q j) exp {(fs (q) + (^4 (k)) + {k, q) exp {-i {y - x) {k - q)) x 

X exp (vJ3 {k) + If4 (q)) exp (f2 (k) - f2 (g)) exp (^23 (g) + f4 (fc))} | 0 ) + ....= 

= - J dk J dg^^^^^^exp(-i(?/-a:)(fc-g)) + ... (12) 


We see that to leading order in the small density ex¬ 
pansion we have that {p {x) p{y)) = p^ — (h^ (x )^(y))i 


[3J|. For example this means that at zero temperature: 
(p(x)p(y)}=p^-^i^^i^^ + .... 

Probability of forming an empty interval. Of particular 
interest interest is the expectation value of the operator 
(exp (aQxy)) in the limit when a —>■ — 00 . In that case the 
only terms that contribute to the expectation value are 
those when there are no particles in the interval [x, y\ and 
Qxy = 0. We call this expectation value the probability 
of having an empty interval. In this limit the expression 
given in Eq. (Uni) greatly simplihes d: 


Pix,y) 


({/c}|exp iaQxy)^^_^ |{fc}) 

m I {k}) 

(0|det(l + ^MT) |0) 
det (1 - ^Kt) 


(13) 

(14) 


Here 


Mt (k, g) = -cf (g) 


exp {^k{y -x) + ^</>(fc)) 
{k — q){k — q + ic) 


X exp ^-^g {y-x)- (g) ) -f 

exp {^q{y-x) P \(t){q)) 

(g — A:) (g — fc -I- ic) 

( i 1 

xexp(--fc(y-a;)- -^(A:) 


(15) 


Here (j)[k) = P (k) + Q {k) with P (k) |0) = (0| Q (k) = 
0 = [P{k),P{q)] = [Q{k),Q{q)] and [P (A:), Q (g)] = 

( (k-qf-^+c ^ ) • leading order expression for this cor¬ 
relation function (keeping only terms with 2x2 matrices 
or less and noticing that P {x = y) = 1) is given by: 


J 


P{x,y) = l- J iy-x) + - / [y - x) 


2tt 2tt 

1 f f (k) f (g) J {k - g)^ + f exp (z {y - x) {q - k)) - 1 exp (z {y - x) {k - q)) - 1 


u 


2tt 27r 


{k-q)" 


{q — k + ic)^ 


{k- q + icy 


For the correlation functions with respect to the ground state, / {k) = 9 {—kp^ kp], and we obtain that: 


kp 


kl 


r-l 2-1 




(k - g)^ 


/ 


exp (z [kp (y - x)] (g - A:)) - 1 exp (z [kp {y - x)] {k - g)) - 1 






For small (z/ - x) this series gives P {x, y) = 1--^ (z/ - x) + ^ {y - x)"^ { (t ~ + (ij:) + ^)) |- 


General correlation functions. We would like to extend our results to general multipoint correlation functions. 
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The most general correlation function can be written as: 

n—1 

O = 00.±1 (xj) exp (16) 

3=0 

Here 0o (xj) = 1, 0i (xj) = (xj), 0_i (xj) = h{xj) 

and = \\vcLe\.Q {y)h{y)dy. Further¬ 

more let us denote by ao the set of Xj where we have 
a 00 {xj)\ o-i the set of Xj where we have (xj) and a_i 
the set of Xj where we have 0_i (xj). Let the cardinal¬ 
ity of a±x be equal to m. In the supplementary online 
information we show that: 


erator O we have the result that is given in Eq. ([5|). We 
have presented a useful expansion for correlation func¬ 
tions of the Lieb liniger gas at short distances and low 
densities. This result can be used in the future to com¬ 
pute correlation functions for quench problems. We hy¬ 
pothesis that a similar relation is valid for arbitrary in- 
tegrable models. In these models there there are many 
particle types. We conjecture the following formula for 
the correlation functions 


{{k}\o\{k})=Y,- 

n\ 

n 


dki 

27r 


dkn 

27r 


X 


({fc}|0|{fc}) ^ (Q| al-oL |0) 

({fc} I {k}) det (1 - 

(17) 

Where UT(k,q) ^ f (k). Arguing again exactly in the 
same way as following Eq. © we can obtain a series 
expansion for a general correlation function thereby pro¬ 
viding efficient methods to evaluate them a deriving the 
LeClair mussardo formula. 

Conclusions. Here we have presented a LeClair Mus¬ 
sardo like formula for generic correlations of the Lieb- 
Liniger gas. We have proven that for a generic local op- 


X ^ (n (fc,) I (kx, ..kF 

Pi Vf=i / 

(18) 

Here is a sum over particle types. It is of impor¬ 
tance to derive this relation for correlation functions for 
a general integrable models. The authors are also cur¬ 
rently working on extending the formula to multi time 
correlation functions. 
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Supplementary online 
information 

I. LECLAIR MUSSARDO FORMULA 


We begin by reviewing some properties of relativisti- 
cally invariant integrable one dimensional systems (for 
simplicity we focus on the single species of particle 
case- generalizations to multiple species of particles is 
straightforward). Eigenstates of integrable models are 
parametrized by sets of rapidities |{0i}). For finite sized 
systems with periodic boundary conditions these rapidi¬ 
ties satisfy bethe ansatz equations 

e^Pi>)^)^Y[S(.Oj-0k) = l (19) 

Here L is the system size, p{9) = m sinh (6) is the mo¬ 
mentum of the particle (here m is the mass of the par¬ 
ticle) and S {9j — 9k) is the scattering matrix between 
particles j and k. In integrable models the two parti¬ 
cle scattering matrix determines all multiparticle scatter¬ 
ing. In the thermodynamic limit when both the particle 
number and the system size is large it is possible to in¬ 
troduce quasiparticle densities. To do so let us denote 
by Lpp (9) d9 as the number of particles in the interval 
\9, 9 -\- d9], Lph (9) d9 as the number of holes in the inter¬ 
val [9,9 -I- d9] and Lpt (9) d9 as the number of states in 
the interval [9,9 d9] so that pt (9) = pp {9)+ph (9). It is 

also convenient to introduce / (0) = ■ The quasipar¬ 

ticle density satisfies the so called thermodynamic Bethe 
ansatz equations: 

1 rJO' 

Pt is)=^p' (^)+y ^ - ^') Pp (^') (20) 

Here p{9) = —i-^\og{S (9)). This is the continuum 
version of Eq. di). Using these quantities various ther¬ 
modynamic quantities like the Yang-Yang entropy may 
be defined. The Yang-Yang entropy associated with the 
densities, {pp (9 ), p^ (9)}, measures the number of states 


[{/c}) consistent with the densities. It is given by: 


S({P^}) = 

= fZo dO [Ph [e) In (^) + Pp {9) In (^)) ■ 


With these definitions the LeClair Mussardo formula for 
the expectation of an operator with respect to an eigen¬ 
state of the theory may be written as 




dO]^ 


X I (22) 

U=1 


Here f {9) = corresponds to the state |{0}) and 

c (^1) --^n) is the so called connected correlation func¬ 
tions for the operator O. Its definition is given below. To 
define the connected correlator recall the Kinematic pole 
axiom for correlation functions of integrable field theo¬ 
ries 0 (on the infinite interval) which in the simplest 
case says: 

( 0 - 0 ) {9,9[,...9'j, \0\9,9 i,..9n) (23) 

= ( 0 ^( 0 ;-^) -05(0-00) X 
X (^ 9 ^, ...0JY I O I 01, ..0Ar) 


The pole when 0 —^ 0 is universal but the residue 
may be more complex then shown in Eq. ([23]) above. 
As a result of this axiom, the expectation value of 
(01 -I- ei, ...0„ -b £„[ O |0i...0„) has poles when a J, 0. The 
connected part of a correlation function is defined by 
52°,c(^l,-^n) = 


Finite Part 


(^lim (01 -b ei, ...0„ -b e„| O |0i.. 

\ei40 



(24) 


In particular it does not have any factors of the form 



II. LIEB-LINIGER GAS AS A 
NON-RELATIVISTIG LIMIT OF THE 
SINH-GORDON MODEL 


The LeClair Mussardo formula has been defined for 
relativistic field theories. We would like to explain why 
it should work for the Lieb-Liniger gas. We would like to 
note that it has never been rigorously proved for general 
correlation functions. We closely follow the derivations 
given in . In the work the authors show that the Lieb- 
Liniger gas (given by the Hamiltonian in the main text) 
can be obtained as the nonrelativistic limit of the Sinh 
Gordon model. The Sinh Gordon model consists of a sin¬ 
gle massive boson governed by the following Lagrangian: 




(s) (<=»!> W)-l) (25) 
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Here v is the velocity of light. The particles in this 
model have mass 


, 1 sin (tt a) 

4 na 


(26) 


Here a = j • The non-relativistic limit is obtained 

STT + Cp-^ 

by the following procedure 


V —>■ oo, g ^ 0, vg = const (27) 


In this case it is possible to obtain, by ignoring some 
terms oscillating at a frequency the Lieb-Liniger 

Hamiltonian given in the main text with the coupling 
constant 


c = 


2 2 
V 9 


16 


(28) 


We note that the results in this section are not a proof 
of the LeClair Mussardo formula which is merely a hy¬ 
pothesis for the Sinh Gordon model, but are merely mo¬ 
tivational for our discussion below. 


Here |0_,) is the state with no bosons on the interval 
Next we use the identity [^: 


b {Xj) n B (fcf) | 0 j) = -i^ aj {k^j) 

X n n 


i=i 


fef/fcs. 


fef/fcs. 


(Ojl n ^3 H ^3 (ko.j) 

X n n (^1^3 (kf) 


1=0 




(34) 


Next we use the relationship 


i: 


III. MULTIPOINT CORRELATION FUNCTION 


We would like to extend our results to general corre¬ 
lation functions. The most general correlation function 
can be written as: 

n—1 

O = (j)o,±i (xj) exp (djQxj,xj+i) (29) 

j=0 

Here i^o (xj) = 1, (t>i (xj) = (xj), (j)-! (xj) = b{xj) 

and Qxj,xj+i = lim^o b^y)b{y) dy. Because of 

translation invariance we may as well assume xq = 0. 
Furthermore let us denote by ao the set of Xj where we 
have a (j)o (xj)', ai the set of xj where we have b^ (xj) and 
a_i the set of xj where we have 4>-i (xj)- Let the cardi¬ 
nality of a±i be equal to m. Now we want to calculate 

(0| C (fcf) ...C {k%) OB (fcf) ...B (k^) |0) (30) 

Take the limit k^ = k^ and express everything in 
terms of a determinant given in Eq. (HZJ in the main 
text. Now we use the notation: 


n 


T(o,l) (fc) (0,L) = nT’G„.,+o(^) 

j=o 



Aj (fc) 
C, (fc) 


Bj (k) \ 

D, (fc) ) 


(31) 

(32) 


Here Xn+i = L and 7)^ (k) is the transfer matrix for 

the interval (a, b). We now use the relationship 0: 


n 

B (kf) ...B (fc^) |0) = ^ n (kf) 10,) X 

X n (kf) dj (kk) f {kf,k^) 

{)<j<k<n 

n 

{Q\c{k^)...c{k^^) = Y (o.in^^(^f)x 

{fcB}=uy^^{fcf} J=o 
X n dk{kf)a,{k^)f{k^',kf) 

0<j<k<n 

( 33 ) 


muZoCAkf)UZiBikP)\o,) = 

= n,>fe 9 {kf , kf) g {kf, kf) (0| det„^. {S^) |0) 


(35) 


Here 


BL = i {kf, kf) a, {kf) d, {kf) exp ($a^ (fcf) -f $£,. {kf)) + 

+ t {kf,kf) aj {kf) dj {kf) exp {<^>a, {kf) + {kf)) 

(36) 


Here 


'^Aj {k) — Qaj {k) + Boj {k) '^Dj — Qoj {k) + Pa^ {k) 

(37) 

The only nonzero commutation relations are given by: 


Puj {k) , Qdi (g)J = Sji In {h (fc, q)) 

Pa^ {k), Qai (g)| = Sji In (fc {q, fc)) 


The relevant expectation values are given by: 


Pa |0) = (0| Qa = 0 (39) 


We note that the state |0) is not related to the state 
with no bosons in any way but is merely an auxiliary 
vacuum used for the purpose of calculating correlation 
functions. Once again 

[$,(fc),$b(g)] = 0 (40) 


With these results we can see that: 
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(0| C (fcf) ...C {k%) OB (fcf) ...B (k^) |0) = ll g {kf, k^) g {k^kf) ^ x 

X X afe {kf) dj {kf) h {kf,kf) x 

{fe'"}=U"=o{fef}U(ci}{fcfo} o<i<'=<" 

X n (^0 ^ (fcfe > kf) X n (0| det (^^■) |0) x 

0<j<fc<n j- 

X Yl y^aj{kfj) h{kfj,kf) X Yl y^dj {kff h{kf,kff 

{“-1} {“1} 


where {Pb\ is the parity of the permutation 
{1,2,3, ...iV} ^ U”^o U {fcf}) and [Pc] is 

the parity of the permutation (1,2,3, ...A^} —>■ 

yj-^,{kfvj{kf}) , and Sl^ = Now intro¬ 
ducing l - Eie{a_i},z<il- We have 

that 

(0| C {kf) ...C {kf) OB {kf) ...B {kf) |0) = 
llj>k 9 {kf,kf) g {kf,kf) (0| det S 

(41) 

To define S we need the following functional fields 


Furthermore we will need the functional fields: 


IpA^ (k) = QAj (k) + PDj (fc) V'D,- = QDj (fc) + PAj {k) 
PAu (k) = TAk {k) + SDk (k) PDk = ro^ (fc) -h SA^ (k) 

(45) 

Here j G o-i, k G ai. The non-zero commutation 
relations are given by: 

[pDj (k) , qoi (<?)] = Sji In {h {k, q)) (46) 

[PAj (k ), qAi (<?)] = Sji In {h {q, k)) 

[sD,- (fc), roi (g)] = Sji In {h {k, q)) 

[sA^ (fc), TAi (<?)] = Sji In {h {q, k)) 


(fc) = Qaj,^ (k) -h PDj^k (k) (k) + Paj^^ (k) 

'^Aj^k (k) = QAj^k (k) + Poj^k fk) (k) + PAj,k (k) 

(42) 


Here 0 < j < k < n. The nonzero commutation rela¬ 
tions are given by: 


The relevant expectation values are given by: 


PDi,k (k), Qdi,,, (<?)] = Su,k),ii,in) ln(/i(fc,g)) 

PA,,Ak),QA,^ (g)] = ln(/i((?, fc)) 

Poj^u (k) , Qdi,^ (q) = 5 (j,k)Xi,m) In {h (k, q)) 

Pa,,^ (k ), Qa,™ (q)] = S(j,k),ii,m) In {h (q, k)) 


(43) 


The relevant expectation values are given by: 


Pa|0) = (0|(7a=Sa|0) = (0|ra = 0 (47) 

In terms of these variables 

m 

S=J2 P'YfldiA (48) 

j&ao j&ai j&a-l 1=1 


Pa |0) = (0| Qa = Pa |0) = (0| Q, = 0 (44) Here: 


J 


=-A Ail 11®° n^° n^° 11“° ^ 

o>j o<j o>j o<j 

X exp 4 >a,,„ {kf) + Y, {kf) + Y {kf) + ^ Ta„,, {kf) 


^o>j 


o<j 


o>j 


o<j 


Furthermore: 


'S'fcj = Sff X exp (V'A, {kf)) 

SlY = ^ ®xp {pD, {kf)) 


(49) 


( 50 ) 
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We now note that: 

= c Oo (fcf) n do {kf) do (k^) Oo (fcf) X dj (fcf) a, {kf) x exp (fc^) + ipA, (kf)) 


o>l 


o<.i 


0>2 


o<3 


exp Y. H (fcf) + 5 ] + J2 '^Aoj (kp) 


(51) 


o>j 


o<j 


Where j is the d’th entry of ai and i is the d’th entry of a_i. Next we need to manipulate these formulas a little 
bit, we start with: 

N N 

(0| det S' |0) = a (kf) d {kf) (0| exp ($a„ {kf) + {kf)) exp {kf) + „ {kf)) det S |0) 

3^1 j=l j<n 

N N 

= n « (kf) d {kf) n h {kf, kf) (6| det S |0) (52) 

j=l j,k=l 


Where 


N 


(6| = (0| exp (Pz,„ {kf) + Pa,, {kf)) exp ($i,. „ {kf) + 'Saj,,, {kf)) 


(53) 


i=i 


j<n 


and Si^m = Si^md {kf) a {kf) exp {—^A„ {kf) — {kf)) ■ Furthermore (O | O) = 1. Now introducing = 

- (6| |0) , = $£,„ - (6| |0) , - (6| „ |0), - (6| |0) and 

specializing to the case when kf = kf we obtain that [^: 

~Sin.= ( E + E + E ' + E A') exp {kf) - {kf)) a {kf) d {kf) X (54) 

(55) 


1^1 


yi^Q;o,i^n j€o:i j^cx-i 

X n exp {-<i>D,,„ {kf) - {kf)) +Si„,[l + J2k {kf,kf) ) + s; 


N 


Im 


j<n 


Where the expressions for + S^f + S/’ ^ + A' we have „ —>• „ and 4'£) 


Here 


Sfm = {d {kf) a {kf) t {kf, kf) {kf) dn {kf) + a {kf) d {kf) t {kf, kf) a„ {kf) d„ {kf) x 
X exp {kf) + {kf) - {kf) - {kf))) aj {kf) dj {kf) 


(56) 


j<n 


Now we factorize out the part oc Now using the relationship Q: 
dim{L + j:f^,K{kf,kf)). First we note that in 

the thermodynamic limit L + K {kf, kf) = {0\C (^kf) ...C {kf) B {kf) ...B {kf) \0) 

2TTLpt{kf). Now introduce the matrix 0/^ = = O^i 27rApt (fc®) x / (^f > ) det (l - ^ATt) 

5im2T:Lpt {kf). Using this matrix we can see that 


(59) 


(6| det S |0) = det 0im (6| det + -^Ut {k, q) 


| 0 ) 


We obtain that 


Where 




/3i = ..=/3,„=0 


| 0 ) 


Ut {k, q) = f iq) S {k, q) - 6im (^L {kf, fcf 


m I {k}) 


det (1 - ^Kt) 


(60) 

Arguing again exactly as in the main text we can obtain 
(58) a LeClair Mussardo like formula. 
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IV. FIELD FIELD CORRELATION FUNCTIONS 
(EXAMPLES) 


We can consider two cases when /2 < 0 and when p, > 0. 
In the case when /2 < 0 we obtain that: 


A. Field Field correlations 

We would like to give an example of how the vari¬ 
ous formulas to the field field correlation functions work 
in practice. These formulas are useful for calculation 
the short distance correlation functions for various ini¬ 
tial states. For example keeping the leading order term 
in Eq. ([^1) we obtain that: 

(6^ (x) 6 = J dk ^ exp [—i {y — x) k] dk + ... 

(61) 

We would like to evaluate this expression for vari¬ 
ous states / (fc). We will concentrate on thermal states. 
In the limit of zero temperature we have that / (fc) = 
6[—kF-ikp\. From this we obtain that 


(it 


WMri). = E(-'U'vs “P-J 

n—1 \ H Y 


l{x-y) 


n(3 

(64) 


In the case when p > 0 we may approximate k'^ — p = 
9 {k) {k — ^/p) 2^/p — 9 (—fc) (fc -I- 2^/p. Using this 

expression we obtain that: 


(it M i (!,)), = 


(65) 


X 2 F I 


■ \x-y\ 




2 + z- 


x-y\ 

Py/P 


,-e 


-/3/x 


( 66 ) 


where 2 E i is the hypergeometric function. 


(b'f {x)b{y))^ 


1 sin {{x — y) kp) 
TT X — y 


(62) 


B. Velocity probability distribution 


This is a good approximation for short distances and 
low densities. For the case of a finite temperature state 
we have that / (fc) = + O (^) where p = 

/i -I- ^ J dgln ^1 -1- “^)^ Q. From this we obtain 

that 


{b^ {x)b{y))^ = j 


dk exp [—* (y — x) k] 
2 tt exp (/? (fc^ — /2)) -I-1 


+ 0 



The velocity probability density is an easily measurable 
experimentally relevant quantity It is given by: 

Piv)-^ JdxJ dy exp (y - x)^ {b'^ (x) b (y)) (67) 

The factor of ^ in the above formula comes from the 
fact that the mass of the Lieb Liniger gas is normalized 
to i. For the simpler case where f (k) = f (—fc) Eq. ([5]) 
simplifies for both x > y and x < y giving: 


{iP {x)b{y))^ = J dk^^exp [-i{y-x) k]dk + J dk^^ J dq^^"^ exp [-i {y - x) k] + 


Now using the relations: 



dx exp \^—ik |a:| 
|a:| exp ( —ik |a;| 


i^x ) = —P.V. 

.V 

‘2" 


k- 


-P. 




i^x ) = P.U.- 


{k+ir 


-P.V. 




■‘‘n 




( 68 ) 
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We get 




1 f{k) 


2tt 


1 

Yc 


(i) (»(‘' i) (‘ (i' (‘' i)+i) (^' i) - (Y 

(4 (i) (‘'i)) -^ (i) (5'^)+!/-' i)) . <«») 

Where c.c. means the complex conjugate of the entire expression and P.V. is the principle value part. This formula 
becomes very accurate for large v where the principle value part becomes subleading. For this case (if we denote by 
N = J dk^^) we have that 


P,.,),. i;(I) p«_ i,(0,1) (,. (£.0) (o,|) P(0,1)(0,1) P2.- (£.0)) 

i(2,.-.(£,o)-4r-(o,|))}}p2«£t(/(£)r-(o,|))Pc.c. 




